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Abstract 



In this paper the evolution of a quantum system drived by a non-Hermitian Hamil- 
tonian depending on slowly-changing parameters is studied by building an universal 
I high-order adiabatic approximation(HOAA) method with Berry's phase , which is 

■ valid for either the Hermitian or the non-Hermitian cases . This method can be 

regarded as a non-trivial generalization of the HOAA method for closed quantum 
system presented by this author before. In a general situation, the probabilities of 
adiabatic decay and non-adiabatic transitions are explicitly obtained for the evolution 
of the non-Hermitian quantum system. It is also shown that the non-Hermitian ana- 
log of the Berry's phase factor for the non-Hermitian case just enjoys the holonomy 
structure of the dual linear bundle over the parameter manifold. The non-Hermitian 
evolution of the generalized forced harmonic oscillator is discussed as an illustrative 
examples. 



1 accepted for publication in Physica Scripta 



1 



1 .Introduction 

Since Berry's phase factor(BPF) was discovered in evolution of the quantum 
system with adiabatically-changing parameters [1], a few methods studying non- 
adiabatic evolution of a quantum system driven by a Hermitian Hamiltonian de- 
pending on slowly (but not adiabatically)-changing parameters have been presented 
in connection with EPF's [2-5]. In these methods, the high-order adiabatic approx- 
imation (HOAA) method was proposed by this author for the first time [4] and has 
been used and developed for many case [6-12]. However, all of these studies has not 
been concerned with a kind of important quantum system, the quantum open system 
that possesses a non-Hermitian (nH) Hamiltonian. This paper will be devoted to the 
generalization of the HOAA method for such nH quantum system. 

In fact, though the Hamiltonian for a closed quantum system, which is usually 
considered as a basic object in quantum theory, must be a Hermitian operator, many 
theories, such as the Fock-Krylov theorem [13], show the probability to apply the nH 
Hamiltonian for those quantum phenomena with dissipation, decay and relaxation 
[14,15]. Recently, many practical problems including the multiphoton ionization, 
the supermode free-electron laser and the transverse mode propagation in optical 
resonator [16-19] have been concerned with the use of nH Schrodinger equation and 
nH Hamiltonian correspondingly. Since the occurrence of a EPF or its analog may 
be established when something in the considered system is varied, it is natural to 
generalize the concept of EPF for nH quantum system. More recently, some authors 
made this generalization and applied it to concrete physical problems [20-21], but 
their studies were only focused on the adiabatic case that the parameters change 
so slowly that transition between any two (quasi-) energy levels do not happen. In 
this paper we especially emphasized the non-adiabatic evolution of the nH quantum 
system and the geometry of the nH analog of EPF. 

The paper is arranged as follows. In section 2, using a similarity transformation 
of the nH Hamiltonian, we build an universal formalism of the HOAA method for 
the nH quantum system. It is also available to the Hermitian quantum system. In 
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section 3, we explicitly analyse the conditions under which the lowest order approxi- 
mation, namely the the adiabatic approximation, can work well. We also compare our 
results with that obtained with the biorthonormal state method [20-21] in adiabatic 
case. In section 4, we apply the general result to a simple toy model - the nH forced 
oscillator to show the usefulness of our analysis. In section 5, we show that the nH 
analog of EPF appearing in an adiabatic evolution of nH quantum system are non- 
unitary holonomy group element for the dual bundles over the parameter maniford. 
In the further studies, we will provide more applications of this generalized HOAA 
method to some physical problems. 

2. Generalized HOAA Method for nH Quantum System 

Let us begin by setting some notations. The Hamiltonian of the quantum open 
system we consider as fillows is a nH operator 

H = H(t) = H[R] = H[R(t)] 

that depends on a set of slowly-changing parameters 

R=R(t) : (i?i(t),i?2(t),...,i?L(t)). 

We now assume that H(t) is diagonalizable at each instant t, i.e., there exists a 
similarity transformation 

U{t) = U[R] = U[R{t)] 

such that 

ei{t) 

€2{t) 



U{t)}i{t)U{t)-' = Hait) 











eN{t) 



H,(t) (2.1) 



where the "quasi-energy levels" 

efe(t) = efe[R], A; = l,2...,iV 
may be complex and U(t) is not unitary correspondingly. 



Let 

>= u{t)mt) > 

be a solution of the nH Schrodinger equation 

tn^^\^{t)>=U{t)\^{t)>. 2.2) 

The equivalent wavefunction \^{t) > must satisfy an equivalent Schrodinger equation 
(ESE) 

zn^|$(t)>=He(t)|$(0>. 2.3) 
with the equivelent Hamiltonian 

ll^{t)^U,{t)-ihU{tr'^, (2.4) 



Separate He(t) into two parts, the diagonal part 

d 

Uo{t) ^ Ud{t) + diagonal part of [-ihU{ty^ — U{t)] 
and the off-diagonal part 

d 

V(t)= of f diagonal part of [-ihU{t)'^—U{t)]. 

Then, 

He(t)=Ho(t) + V(t). 

Late on we will show that the diagonal part Ho(t) governs an adiabatic evolution of 
the nH quantum system while the off-diagonal part V(t) governs the non-adiabatic 
transitions among the quasi-energy levels. In fact, since V(t) completely vanishes 
when H is independent of time t and H = H(t) = H[R] is a smooth function of R(t) 
, we can regard V(t) as a pertubation when H or R depends on time "slightly", i.e., 
R(t) changes as t slowly enough. 

Based on the above decomposation of He(t), we can use the standard time- 
dependent pertubation theory to solve the ESE (2.3), obtaining a approximate series 
solution 

\^{t) |$0(i) > +|$^(t) > +|$2(t) > +... + \^\t) > +..., (2 - 5) 
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Here, the I'th order solution > is determined by 



d 



dt 



(2.6) 



/ = 1,2,..., 



Because the above equations of Z'th order solution > only concern {I — l)'th 

order solutions > , we can get each order approximate solution starting from 

the zero'th order one 

in Jo 

|$(0) >= ?7(0)-^|^(0) > . (2.7) 
In order to write out thse solutions in a quite explicit form, we use the eigenstates 



|1 >= 
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of Ho(t) with the corresponding eigenvalue 
where the additional term to the quasi-energy 



7n(^) = ^/^ <n\Uitr'^^\n>dt' 



10 dt 

will be proved to be the nH analog of Berry's phase. Then, each order solution 



(2.8) 



N 



|$'(t) >= ^Ci(t)e*T'"We*/o^"(*')'^*'|n >, 

n=l 

follows from eq.(2.6) immediately. Here, the coefficients C^(t) satisfy 

C^{t) =< n|$(0) >=< n|C/(0)-^|*(0) >, 



(2.9) 



Cl{t)^Y.-^ l <nm')\ 

m=l •'^ 




c'-\t')dt', 



(2.10) 



where 



n,m — 



^n{t) - 6m{t) + 7m (0 - 7n(^) 

h 



Notice that the difficulty that the eigenstates of H are not orthogonal to each 
other due to non-Hermiticity of H has been avoided in the above discussion by a 
trick building the ESE (2.3) to find the pertubation V(t). This is the key to our 
studies in this paper. 

S.Adiabatic Approximation and Comparision with Biort honor mal State 
Method 

In this section our focus first is on the adiabatic conditions that the zero'th order 
(adiabatic) approximation can well approach the true evolution of the nH quantum 
system. Consider the integral 





Jo 



appearing in the first order approximation 




m=l 



< n|V(t')|m > e'Jo 




Here, we have separated the damping factor 




and the oscillating factor 



,iRe[Jl^ u!„,m{s)ds] 



If the latter oscillates so fast that the conditions 



<n|V(t)|m>e-''"[/o 
Re[u;n,m{t)]h 



is)ds] 



< 1 



(3.1) 
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hold for m ^ the integral I tends to zero. This statment can be proved by 
integrating I by part. Then, we get so-called adiabatic conditions (3.1). When they 
hold for a nH quantum system, we can ignore the higher- approximation solutions 
\<^i{t) >,(Z = 1,2,...) . 

To compare our analysis with the biorthonormal state method in refs. [20,21] in 
the adiabatic case, we define 

\^n{t) >=|0„[R]>= U{t)\n>, 
\Xn{t) >= |xn[R] >= [U{t)-']^\n >, (3.2) 

Obiviously, they are the eigenstates of H(i) and H(t)^ respectively with the eigenval- 
ues e„(t) and e*(t). Using the completeness relation and the orthonormal relations 
of states \n> {n = 1, 2, A^): 

N 

\n >< n\ — I{unit operator), < n\m >— 5m,n, 

n=l 

we immmediately get the generalized completeness relations 

N N 

E \Xn{t) >< Mt)\ = E \Mt) >< Xn{t)\ = I (3.3) 
n=l n=l 

and the biorthonomarl relations 

< (f>m{t)Xn{t) > = < XmitMnit) >= ^m,n- (3.4) 

In terms of the biorthonormal basis {|^„(t) >, |x„(t) >}, we reexpress the high- 
order approximations of the nH quantum evolution as 

oo oo 

m) >=Ei*'w >^j:u{tmt)> 

1=0 1=0 

oo N ^1 

= E E Cl{t)e^^-^'^e^nh -(*')'^*'|0,(t) >, (3.5) 

1=0 n=l 

with the coefficient equations: 

C°(t) =<Xn(0)|^(0)>, 

cut) = - E ^ / < Xn{t')\Mt') > --'-^^)^^CU\t')dt', (3.50 

m=l 



for / = 1, 2, 3, .... Here, the nH analog of Berry's phase is rewritten as 

7„(t) ^ij^ < Xn{t')\j^(t>n{t') > dt' (3.6) 

Correspondingly, the adiabatic conditions (3.1) are rewritten as 

, <Xn(0l|7c/>.(t)>e-^<-"--(-)'^-] 
Re[u;n,m{t)] 

Now, let us consider the adiabatic evolution of an instantaneous eigenstate of H(t) 
under the adiabatic conditions, if the system is initially in state 

|*(0) >= U{t = 0)\n>=\<j)n{0) >, 

|$(t) >= e*^"We^/o ^"(*')'^*'|0„(t) > 

= e^""W|0„(O> (3.8) 



it will evolve into 



where 



^n{t) = ln{t) - I fen{t')dt'. 

n Jo 



The above eq. (3.8) manifests that the adiabatically- evolving state is always an 
eigenstate of the instantaneous Hamiltonian H(t) if the initial state is such state at 
t=0. This is the quantum adiabatic theorem for nH quamtum systems, which is quite 
similar to that for the Hermitian case. It shows the quantum number labelling the 
quasi-energy level to be an adiabatic invariant. Since 7„(t) and e„(t) are usually not 
real due to the non-Hermiticity of H(t), the damping factor ea:p[— /mQ„(t)] causes 
the adiabatic decay of the state. Especially, when e„(i) is real, the decay only results 
from the nH-analog of BPF and is a purely geometrical effect. In next section an 
example will be used to illustrate this situation. 

When the adiabatic conditions (3.7) are broken, there may appear the transitions 
from an instantaneous eigenstate |0n(O) > to others |0„(t) > for m 7^ n. From 
eqs.(3.5) or (2.5,2.9,2.10), we obtaine the transition probabilities 

P(n ^m) = h t < m\Vit')\n > -^As)ds ^^f^2^-2im[m] (3.9) 
% Jo 



or 



(3.90 



Therefore, we reach the statement mentioned in section. 2 that the diagonal and off- 
diagnoal parts respectively describe the adiabatic evolution and the non-adiabatic 
transitions for the quantum system. Notice that Berry also discussed the geometrical 
phase effect in non-adiabatic transition in connection with the non-Hermiticity due 
to the complexization of time [22] . 

4. Example: Generalized Force Hamornic Oscillator 

Now, let us apply the above general analysis to a simple model - the generalized 
forced harmonic oscillator. The time-dependent Hamiltonian for this model is 



where a — a{t) and P — j3{t) are slowly-changing complex parameters, the constant 
uj is real , a and a) are respectively the creation and annihilation operators for a 
boson state, which satisfy 



For the usual forced harmonic oscillator (FHO), there is constraint a = j3* and the 
corresponding Hamiltonian is Hermitian. However, we make here a generalization to 
remove this constraint so that the Hamiltonian (4.10) is non-Hermitian. 
Using the translated boson operators 



H(t) = H[q;, P] = %uj{a}a + (5a^ + aa) 



(4.1) 



[a, a'] = 1. 



A{a)'^ ^ + a, A{(5)^a + I3 



(4.2) 



obeying the same boson commutation relation as the above 



[A(/3),A(a)t] = l, 



we rewrite down the Hamiltonian (4.1) as 



H(t) = H[q;,/3] = huj[A{a)^ A{I3) - a(3]. 
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Based on the above expression for H(t), we immediately obtain the instantaneous 
eigenstates of H(t) = H[q;,/?] 

\(j)n{t) >= |(/)n[a,/9] >= \n[a,/3] > 

= -^[A{a)rW) > (4.3) 
Vni 

with the eigenvalues 

€„(t) = €„[q!,/3] = (n — a/3)fiu!. 
Here, the translated vacuum state |0(/3) > obeys 

AiPMP) >= 0, or a\QiP) >= -/?|0(/?) >, 

that is to say, |0(/9) > is a coherent state 

W)>^C^,f,e-^'''\0>, (4.4) 

where |0 > is the usual vacuum state satisfying a\0 >= and Ca,i3 is to be determined 
by the biorthonormal conditions . According to the above analysis, we can finally 
write the explicit form of |0„(i) > 

fe=o {n-k)Wk\ 



where 



\k >= ^[a¥|0 > 



is the k boson state. 
Noticing 

H[a,/3]t = H[P*,a*], 
we easily obtain the dual vector |Xn(t) > to |0n(t) >: 

\Xn{t) >= \Xn[a,P] >= HP*, a*] > 

as the eigenvectors of H[q;, /?]''' with the corresponding eigenvalue en[(3*,a*]. Then, 
the biorthonormal conditions 

< Xn{t)\(l)m{t) >= Sm,n 
10 



define the biorthonormalization coeeficient 

Taking into acount the translation transformation 

e^'^a^e-^" = + 

e^"^ae-^"^ = a - (4.5) 
^ G complex field C 
we get the U(t)-operator expression of the biorthonormal basis >, \Xn{t) >} 

\4>n{t) >= |0n[R] >= U{t)\n>, 

\Xn{t) >= |Xn[R] >= [U{t)-^\^\n >, (4.6) 

where the U(t)-operator 

U{t) = /?] = e-5°/3e-^«^e-"" (4.7) 
diagonahzes the Hamiltonian H(t), i.e., 

U{t)-^B.{t)U{t) = fiuja)a. 

Now, let us discuss the evolution governed by the nH Hamiltonian (4.1). To this 
end, we first calculate 

+]^[a(5 - (5a]5m,n (4.8) 
Then, we obtain the nH analog of EPF 

2 Jo 

= ea;p{- ^^^j!^d(3 - (3da\] (4.9) 
where C is a curve {R(t) } in a four dimensional parameter manifold 

M :{R = {a, 13) = (cti, a2, fii, I32)\a = cti + ia2,l3 = /3i + i/32} 
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where ai,a2,(3i and (32 arc real. Because the factor e*^"*^*^ is independent of the 
quantum number n, the adiabatic evolution 

N 

\m >= e'^W J2 < Xn(0)|^(0) > e-'--'\Mt) > 

n=l 

N 

= e-rWe^®W Y: < Xn(0)|*(0) > e-'^'^'lMt) (4-10) 

n=l 

is accompanied with a global geometrical factor with the oscillating part e*®^*'' : 

e{t) = Re[-fn{t)] = J /*[/3id2 + p2ai - a2$i - ai$2\{t')dt' 
2 Jo 

and the damping part e~'"^*^: 

T{t) = Im[jn{t)] = ^ I\(32a2 - A«i - "2/32 + ai(3i]{t')dt' 

Z Jo 

Notice that, when a{t) = p{t)* , the BPF of coherent state for FHO in ref.[23] is 
given once again as a special case of eq.(4.9). 

For the case with large quantum number n or rapidly-changing parameters a{t) 
and /3{t), the adiabatic conditions 

Vnm\e-m « 1- V^\m^-m « 1 (4.11) 

do not hold and we need to consider the non-adiabatic effects caused by the first 
order approximation at least. If the system is initially in a state |(/>fc(0) >, the initial 
conditions 

C°(0) = 5„,, , C»=0,/or l>l 
leads to the first order approximate solution 

Jo 

+^fcne-^^We*[®W-(*^+^)'^*] t Kty^UtUk+iit) > (4.12) 

^0 

The probabilities of the transition from |0„(O) > to |0„(t) > (n=k-l, k+l) are 

Pik ^ n) = fce-2rW| f a(t')e-'^'' dt' \Hnk-i + Ije-^^Wl [' p(t')e'^''Ut\^6nk+l■ 
Jo ' Jo ' 

(4.13) 
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Obviously, the selection rule for such transition is 

An = +1 Or - 1, (4.14) 
under the first order approximation. 

5. The nH- Analog of BPF and Holonomy in Dual Line Fiber Bundles 

For Hermitian quantum system, Simon recognized that the BPF is precisely the 
holonomy in a Hermitian line bundle defined by the adiabatic evolution. The adia- 
batic evolution can be interpreted as a parallel translation in such bundle [24]. Now, 
a question naturally arises for the non- Hermitian case: what is a geometrical inter- 
pretation of the nH analog of BPF? The answer will be given in this section. 

Let M be the parameter manifold formed by the parameters R = (i?i, R2, Rl) 
. A line bundle defined by the nH Hamiltonian H[R] is 

F„ = {(R, |(7„[R] >) I H[R]|(7„[R] >= e„[R]|(7„[R] >,R e M}. 

Its fiber space is a one-dimensional linear space 

Vn : {|(T„[R] >= e*^'^^'0„[R] > | ^[R] is a real functions depending on R} 

Over the same base manifold M, the dual bundel F* is defined by 

K = {(R, >) I H[R]t|<[R] e:[R]|<[R] >,R e M}. 

where the fiber space 

V* : {|cr*[R] >= e'^'[^l|x„[R] > | ^'[R] is a real functions depending on R} 
is the dual space of Vn and |Xn[R] > is the dual elemet of |0n[R] > , i-e., 

< Xn[R]|0n[R] >= 5m,n- 

Since the quantum number n labelling an eigenstate |0n[R] > of the instantaneous 
Hamiltonian H[R] is an adiabatic invariant, we can assume that 

Wnit) |(7„[R] >= C„[R]|0„[R] > 

13 



is an evolution state in adiabatic case. Now, let us show that the holonomy group 
elements of F„ ( F*) is the nH analog of BPF for the adiabatic evolution while 
{Kit) > \t e [0,T]} ( {\a *n {t) > \t e [0,T]} ), as a curve in F„ (F;) , is a 
horizontal hft of the curve C : {R(t)|i e [0, T]} on base manifold M. To this end we 
consider a decomposition of a tangent vector 

ft\^n{t) >^ f^Kit) > +f^Mt) > (5.1) 

where the vertical part along the fiber is 

^Wn{t) > = < Xn[R]|^kn[R] > |0n[R] > 
= (^a[R]+ < Xn[R|^|0n[R] > C„[R])|0„[R] > (5.2) 

and the horizontal part orthogonal to the fiber is 

^K{t) >= E < Xm[R]4kn[R] > |0m[R] > . (5.3) 



A parallel-translation implies 

dt 



<Xn[R|^k„(i)>=0, (5.4) 



which results in a one-form equation 

rfC4R]+ < Xn[RM<^n[R] > C„[R] = 0. (5.5) 

Its solution 

C„(t) = C„[R(t)] = e^^"W 

just gives the nH analog of BPF again. 

For a cyclic evolution that R(0) = R(t) and C is a closed curve, the complex 
phase 7n(t) in the nH analog of BPF can be expressed as an closed path integration 

7„(r)=7,[C]=£A„[R] (5.6) 

of a complex potential one-form 

A„[R]=Z<Xn[R]M<^n[R]> 
14 



Similarly, on dual bundel F* , the parallell-translation condition 

<<^„[R]|^|a*„ (t) >=0 
results in the dual nH analog of BPF with the geometric phase 

%{t) =^jl< UW)]\^,Xnm'^')] > dt'- (5.7) 

Obviously, the nH analog exp[yn{t)] and its dual exp[%(t)] occuring in the adia- 
batic evolution are holonomy group elements on the line bundel F„ and its dual F* 
respectively. Except for the effect of dynamical factors 

the adiabatic effect in quantum evolution are equivalent to the parallel-translation 
on the hue bundels F„ and its dual F* . This is a circumstance similar to that for 
the Hermitian quantum system. 
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